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^ Abstract 

We present data for the axial coupling constant qa of the nucleon obtained in lattice QCD with 
D , two degenerate flavors of dynamical non-perturbatively improved Wilson quarks. The renormaliza- 



tion is also performed non-perturbatively. For the analysis we give a chiral extrapolation formula 



^ ' for qa based on the small scale expansion scheme of chiral effective field theory for two degener- 

ate quark flavors. Applying this formalism in a finite volume we derive a formula that allows us 
to extrapolate our data simultaneously to the infinite volume and to the chiral limit. Using the 
additional lattice data in finite volume we are able to determine the axial coupling of the nucleon 
in the chiral limit without imposing the known value at the physical point. 
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I. INTRODUCTION 



The axial coupling constant gA of the nucleon has been studied theoretically as well as 
experimentally for many years. It can be defined as the value of the axial form factor of 
the nucleon at vanishing 4-momentum transfer. More explicitly, one considers the isovector 
axial current A'^~'^ = U'j^'y^u — d'jf^'j^d, where u and d denote the up and down quark fields, 
respectively. We work in the limit of exact isospin invariance, i.e. u and d quarks are assumed 
to be degenerate in mass. The proton matrix element of this current has the form factor 
decomposition 

{p\s'\A;-''\p,s) = u{p',s') 

Here q = p' — p denotes the 4-momentum transfer and u{p, s) is the proton spinor for 
momentum p and spin vector s. The states are normalized according to {p',s'\p,s) = 
{2nY 2p^5{p — p')5ss' , we take = — m^, and mjy is the nucleon mass. So qa = Ga{0) is 
determined by the forward matrix element 

{p,s\Al-''\p,s) = 2gAS,. (2) 

In parton model language, the forward matrix elements of the axial current are related to 
the fraction of the nucleon spin carried by the quarks. Denoting by Am (Ad) the contribution 
of the u (d) quarks, one has 

(p, s|m7^75m|p, s) = 2Aus^ (3) 

and similarly for the d quarks. Thus we can write qa = Am — Ad. 

In this paper we report on new results for qa obtained by means of Monte Carlo simula- 
tions of lattice QCD with two dynamical quark flavors. While Eq. ([2]) lends itself immediately 
to an evaluation on the lattice, it is not yet possible to perform the simulations at physical 
quark masses. Moreover, the size of the lattice is necessarily finite and in practice not very 
large. Thus, apart from the unavoidable continuum extrapolation, we have to cope with the 
thermodynamic limit and the extrapolation to small (physical) quark masses, the so-called 
chiral extrapolation. Guidance for these extrapolations is provided by chiral effective field 
theory (ChEFT). 

In its standard form, ChEFT describes low-energy QCD by means of an effective field 
theory based on pion, nucleon, . . . fields taking into account the constraints imposed by 
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(spontaneously broken) chiral symmetry. As long as one stays in the p-regime of CliEFT 
(with appropriate boundary conditions), the Lagrangian does not depend on the volume. 
So besides the quark-mass dependence the very same Lagrangian governs also the volume 
dependence, and finite size effects can be calculated by evaluating the theory in a finite 
(spatial) volume. Thus the finite volume does not introduce any new parameters and the 
study of the finite size effects yields an additional handle on the coupling constants of 
ChEFT. 

There are several ways to treat baryons in ChEFT. Here we apply the (non-relativistic) 
small scale expansion (SSE) which uses explicit pion, nucleon and A(1232) degrees of 
freedom, and extend the previous calculations of the quark-mass dependence of gA jsl, S| in 
this scheme to finite volume. The dependence on the lattice spacing a could be included 
(for a review see Ref. j^), but we shall not consider this possibility here. 

A somewhat more phenomenological approach to the chiral extrapolation of qa (and other 



nucleon matrix elements) has been developed in Ref. I™. The volume dependence of such 

n n Q 

matrix elements has also been studied by several methods 171, ISI, 19| and we shall compare our 
procedure in some detail with that of Beane and Savage [7| . 



Preliminary results of our investigation have been presented in Refs. 



10 



II. THE SIMULATIONS 



The QCDSF and UKQCD collaborations have generated ensembles of gauge field con- 
figurations using Nf = 2 non-perturbatively 0{a) improved Wilson quarks and Wilson's 
plaquette action for the gauge fields. The simulation parameters are listed in Table [T] along 
with some auxiliary results needed later on. Note that we have two groups of three ensem- 
bles each which differ only in the volume (simulations 9, 10, 11 and 12, 13, 14). In this 
paper we shall not consider any partially quenched results, hence we set k = Hsca- 

When we want to compare (fit) our results with formulae from ChEFT we need all 
numbers in physical units, i.e. we have to fix the scale. This is usually done with the help 
of the force scale [l^ derived from the heavy-quark potential. From phenomenology 
we know that tq ~ 0.5 fm in the real world. Unfortunately, it is not easy to determine 
this value precisely because it is related to measurable quantities only through the use of 
potential models. 
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TABLE I: Simulation parameters together with results for the force scale rg and the pion mass m. 
in lattice units. Also the value of the plaquette Uq is given. 





Coll. 


/3 




'^sea 


volume 


ro/a 








1 


QCDSF 


5.20 


0. 


,1342 


163 


X 


32 


4.077(70) 


0.5847(12) 


0, 


,528994(58) 


2 


UKQCD 


5.20 


0. 


1350 


163 


X 


32 


4.754(45) 


0.4148(13) 


0, 


.533670(40) 


3 


UKQCD 


5.20 


0. 


1355 


163 


X 


32 


5.041(53) 


0.2907(15) 


0, 


.536250(30) 


4 


QCDSF 


5.25 


0. 


1346 


163 


X 


32 


4.737(50) 


0.4932(10) 


0, 


.538770(41) 


5 


UKQCD 


5.25 


0. 


1352 


163 


X 


32 


5.138(55) 


0.3821(13) 


0, 


.541150(30) 


6 


QCDSF 


5.25 


0. 


13575 


243 


X 


48 


5.532(40) 


0.2556(5) 


0, 


.543135(15) 


7 


UKQCD 


5.29 


0. 


1340 




X 


32 


4.813(82) 


0.5767(11) 


0, 


.542400(50) 


8 


QCDSF 


5.29 


0. 


1350 


163 


X 


32 


5.227(75) 


0.4206(9) 


0, 


.545520(29) 


9 


QCDSF 


5.29 


0. 


1355 


123 


X 


32 










10 


QCDSF 


5.29 


0. 


1355 


163 


X 


32 


- 


- 




- 


11 


QCDSF 


5.29 


0. 


1355 


243 


X 


48 


5.566(64) 


0.3269(7) 


0, 


,547094(23) 


12 


QCDSF 


5.29 


0. 


1359 


123 


X 


32 










1 O 

16 


QCDSl* 


5.29 


0. 


1359 


163 


X 


32 










14 


QCDSF 


5.29 


0. 


1359 


243 


X 


48 


5.88(10) 


0.2392(9) 


0, 


,548286(57) 


15 


QCDSF 


5.40 


0. 


1350 


243 


X 


48 


6.092(67) 


0.4030(4) 


0, 


,559000(19) 


16 


QCDSF 


5.40 


0. 


1356 


243 


X 


48 


6.381(53) 


0.3123(7) 


0, 


,560246(10) 


17 


QCDSF 


5.40 


0. 


1361 


243 


X 


48 


6.714(64) 


0.2208(7) 


0, 


,561281(8) 



Yet this is not the only problem with setting the scale. In the analysis of quenched 
simulations it is common practice to identify tq as extracted from the Monte Carlo data 
with the physical tq. However, the latter refers to a world where the sea quarks have their 
physical masses while in the quenched simulations they are infinitely heavy. Still, for many 
quantities quenched calculations lead to results that agree surprisingly well with experiment. 

In simulations with dynamical quarks additional difficulties arise. In this case ro/a de- 
pends not only on j3 as in the quenched approximation, but also on n. Then one can either 
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TABLE II: Bare and renormalized results for g^. The values for L and in physical units have 
been calculated using ro = 0.467 fm together with ro/a at the respective quark masses. 





13 






volume 


L[fm] 


m^[GeV] 


„barc 


9A 


1 


5.20 





1342 


16'^ 


X 32 


1.84 


1.007(17) 


1.452(17) 


1.185(16) 


2 


5.20 





1350 


16^ 


X 32 


1.57 


0.8332(83) 


1.514(29) 


1.201(24) 


•i 

KJ 


5 90 


n 




16^ 


X 32 


1 48 


n 61 Q9l'73'l 


1 3Q6l'36'l 


1 n88('99'l 


4 


5.25 





1346 


16^ 


X 32 


1.58 


0.987(11) 


1.442(13) 


1.176(12) 


5 


5.25 





1352 


16^ 


X 32 


1.45 


0.8295(93) 


1.438(20) 


1.148(17) 


6 


5.25 





13575 


243 


X 48 


2.03 


0.5975(45) 


1.456(10) 


1.1398(98) 


7 


5.29 





1340 


16^ 


X 32 


1.55 


1.173(20) 


1.437(12) 


1.207(12) 


8 


5.29 





1350 


16^ 


X 32 


1.43 


0.929(13) 


1.409(12) 


1.143(11) 


q 


5 29 


n 


1 355 


12^ 


X 32 


1 01 




1 1 81 (PiO) 


n QdyfAR) 


10 


5.29 





1355 


16^ 


X 32 


1.34 


_ 


1.364(25) 


1.087(21) 


11 


5.29 





1355 


243 


X 48 


2.01 


0.7688(90) 


1.459(11) 


1.163(11) 


12 


5.29 





1359 


12^ 


X 32 


0.95 


- 


0.97(10) 


0.763(79) 


13 


5.29 





1359 


16^ 


X 32 


1.27 




1.253(45) 


0.985(36) 


14 


5.29 





1359 


243 


X 48 


1.91 


0.594(10) 


1.413(22) 


1.111(18) 


15 


5.40 





1350 


243 


X 48 


1.84 


1.037(11) 


1.4737(74) 


1.2234(88) 


16 


5.40 





1356 


243 


X 48 


1.76 


0.8420(72) 


1.451(11) 


1.180(11) 


17 


5.40 





1361 


243 


X 48 


1.67 


0.6264(63) 


1.410(20) 


1.127(17) 



identify the value for tq obtained at the given [jS, k) combination with the physical tq or one 
can attempt an extrapolation of r^/a to the chiral limit, e.g. at fixed /3, and use the extrapo- 
lated value to set the scale. The first procedure is analogous to the standard approach in the 
quenched case and ensures a continuous connection with the quenched results as the quark 
mass tends to infinity. The second procedure, on the other hand, has the virtue of leading 
to a mass-independent lattice spacing a. It is however to be noted that, strictly speaking, 
the chiral extrapolation should be performed at fixed = (1 -|- hgam)gQ (see Ref. [l3|) and 
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not at fixed (3, although the difference seems to be rather small in practice. 

We have determined tq in physical units from the nucleon mass. Of course, this can 
only be done after a chiral extrapolation to the physical pion mass. Encouraged by our 
comparison of nucleon mass data in different volumes with chiral perturbation theory fl^ . 
we have used the extrapolation procedure described in [l^ as Fit 1 to recent nucleon masses 
obtained by the CP-PACS and JLQCD collaborations along with updated masses from the 
QCDSF-UKQCD collaboration. Varying the assumed physical value of ro one can make the 
fit curve pass through the physical point, which happens for tq = 0.467 fm. This number is 



consistent with a recent lattice calculation of /^r [15|, giving ro = 0.475(25) fm. A similar 
result for ro was also quoted in [l^ taking as input level splittings in the T spectrum. In the 
following we shall use ro = 0.467 fm, but for comparison we shall also consider ro = 0.5 fm. 

We compute qa from forward proton matrix elements ([2]) of the flavor- nonsinglet axial 
vector current at p = 0. The required bare matrix elements are extracted from ratios of 
3-point functions over 2-point functions in the standard fashion. Our results for g^'^'^ are 
collected in Table [Tll Compared to the computation of hadron masses, additional difficulties 
may arise in the calculation of nucleon matrix elements such as qa'- In general there are 
quark-line disconnected contributions, which are hard to evaluate, the operators must be 
improved and renormalized etc. Fortunately, in the limit of exact isospin invariance, which 
is taken in our simulations, all disconnected contributions cancel in qa-, because it is a flavor- 
nonsinglet quantity. So we do not have to worry about any disconnected contributions. 

The axial vector current is 0{a) improved by adding the derivative of the pseudoscalar 
density with the appropriate coefficient: 

^r^(a^) = ^(a;)7/i.75g(a;) + ac'^^d^q{x)^^q{x) . (4) 

Hence the improvement term does not contribute in forward matrix elements such as ([2]). 
We have denoted the improvement coefficient (usually called ca) by in order to avoid 
confusion with a coupling constant appearing later. 

The renormalized improved axial vector current can be written as 

A, = ZA{l + hAam)^^^ (5) 

with the bare quark mass m = (I/k — 1/kc) /(2a). The values of used here can be found 
in Table IIIII The coefficient hA is required to maintain 0(a) improvement also for non- 
vanishing quark mass. We are not aware of a non-perturbative evaluation of this coefficient 
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for our action. Hence we have to resort to perturbation theory. A one-loop calculation 



yields 



17t 

bA 



1 + b^j^g', + 0{g',) , 6^ = 0.11414(4) Cf , (6) 



where Cp = A/3 in QCD and qq denotes the bare coup 



tadpole improved version of this result, i.e. we take 18| 



ing with P = G/qq. We shall use the 



(l H- i) I) . (7) 

Here uq is the fourth root of the expectation value of the plaquette with the perturbative 
expansion 

uo = l- ^gl + 0{g',) . (8) 

While the coefficient bA will be computed in tadpole improved one-loop perturbation 
theory, we calculate the renormalization factor Za non-perturbatively by means of the Rome- 
Southampton method jiol. I20I1. Thus Za is first obtained in the so-called RI'-MOM scheme. 



Using continuum perturbation [21[ theory we switch to the MS scheme. For sufficiently large 
renormalization scales fi, Z^ should then be independent of fi. However, unless /i ^ 1/a 
lattice artefacts may spoil this behavior. Since our scales do not always satisfy this criterion, 
we try to correct for this mismatch by subtracting the lattice artefacts perturbatively with 
the help of boosted one-loop lattice perturbation theory. As Fig. [1] shows, some lattice 
artefacts still remain, but we can nevertheless estimate Za- As our central value we take the 
result obtained at /i^ = 5 GeV^ while the difference to the result at = 10 GeV^ gives us an 
estimate of the (systematic) error. Adding the statistical error in quadrature we arrive at the 
numbers given in Table UTTt where we also compare our results with a recent determination of 



Za by the ALPHA collaboration [2^. While the numbers differ significantly at the lower (3 
values there is a tendency towards agreement as (3 grows. This suggests that the differences 
are mainly due to lattice artefacts and will not influence the continuum limit of renormalized 
quantities. 

Our results for g^ renormalized in the way just described are given in the last column of 
Table [IT] and are plotted in Fig. [2l Note that the error on gj^ is dominated by the statistical 
error of the bare values while the uncertainty of Za contributes only comparatively little. 
We have taken m^r from the largest available lattice at each (/5, k) combination. The scale 
has been set by means of the force parameter tq with tq = 0.467 fm, and ro/a has been taken 
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TABLE III: Values of Zj^ from this work and from the ALPHA collaboration. The latter numbers 
were obtained from Eq. (3.6) in Ref. [24], where an error decreasing from 0.01 at /? = 5.2 to 0.005 



at /3 = 5.7 is ascribed to them. In addition we give and the values of ro/a extrapolated to the 



chiral limit (from 



2^) 



5.20 



5.25 



5.29 



5.40 



Z\ (this work) 
Za (ALPHA) 



(^o/a; 



ch.I. 



0.765(5) 
0.719 

0.136008(15) 
5.455(96) 



0.769(4) 
0.734 

0.136250(7) 
5.885(79) 



0.772(4) 
0.745 

0.136410(9) 
6.254(99) 



0.783(4) 
0.767 

0.136690(22) 
7.39(26) 



0.84 - 




0.76 




<>00000000000000 0" 

®°°OOOo O CD O O, O O CD 



10 



20 



30 



40 



50 



[GeV^] 



FIG. 1: The renormalization constant Za plotted versus the square of the renormalization scale [i. 



at the given quark mass. Obviously there are considerable finite size effects. A qualitatively 
similar volume dependence has already been observed in quenched simulations {2^. Note 
that the "large volume" results for qa obtained at our smallest quark masses (m^r ~ 600 MeV) 



for all four /5 values lie very close together indicating that discretization effects are sma^ 
Our results may be compared with other evaluations of gA in dynamical simulations 



1. 



261 . l27l |. They all differ in the lattice actions employed: In Ref. [25|] unimproved Wilson 



fermions are used, the RBCK col 
most recent LHPC investigation 



aboration 26|] works with domain wall fermions, and the 
chooses a hybrid approach with domain wall valence 



quarks on improved staggered sea quarks. All these studies obtain a rather weak quark- 
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mass dependence of gA in agreement with our findings. However, most of their values 
he somewhat above ours. At the moment, the reason for these discrepancies is not yet 
clear. One possibility is the different treatment of the renormalization. The perturbative 



renormalization employed in Ref. 



25| entails a relatively large uncertainty because it relies 



on a one-loop calculation only and in general lattice perturbation theory does not converge 
very rapidly. References H] exploit the fact that in the framework of domain wall 
fermions there are five-dimensional (partially) conserved axial and vector currents. One 
should however keep in mind that they are exactly conserved only in the limit of an infinite 
fifth dimension, while the actual simulations are necessarily performed with this extension 
kept finite. 

Note that using the results of the ALPHA collaboration we would have obtained even 
smaller numbers. As indicated above, this difference may be expected to disappear in the 
continuum limit. However, from the comparison of our results at different values of (3 we 
have the impression that lattice artefacts are rather small in the data renormalized with our 
Za values (see Fig. [2]). This is less so when the Za of the ALPHA collaboration is used. 
In any case, there is still some work to be done before we can be sure that all systematic 
uncertainties are under control. 



III. CHIRAL EFFECTIVE FIELD THEORY (INFINITE VOLUME) 

At low energies one can evaluate the matrix element ([1]) of the isovector axial current of the 
nucleon within a low energy effective theory of QCD utilizing the methods of ChEFT. In the 
forward limit (Eq.([2])) one then obtains the axial coupling constant qa- In the SSE formalism 
the results appear as expansions in powers of a small parameter e, which collectively denotes 



small pion four-momenta, the pion mass, baryon three-momenta and the nuc 



eon-A mass 



.putt., . tKe 0... ....t Q. At one Ka,s to evaluate S d,a,a.s Q, wMeK a.e 

displayed in Fig. [3|, involving nucleon, pion and A(1232) degrees of freedom. When ChEFT 
calculations of nucleon properties using just pions and nucleons as the active degrees of 
freedom are extended to ChEFT calculations of the same quantities also employing explicit 
A(1232) degrees of freedom, one obtains renormalizations of the chiral limit couplings of 
the nucleon by polynomial terms ~ (Aq)", where Aq denotes the (finite) N-A mass splitting 
in the chiral limit. This would naively imply that the chiral limit properties of a nucleon 
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FIG. 2: Results for ^a- For (/?,«;) = (5.29,0.1355) and (5.29,0.1359) we show the results obtained 
on three different spatial volumes. The smallest box size L, leading to the smallest qa, is about 
1.0 fm, while the middle size is approximately 1.3 fm. All other volumes are larger. The star 
represents the physical point. 



differ between ChEFT schemes with and without resonance degrees of freedom. In the SSE 
formalism such an unphysical scenario is automatically avoided because Aq (by construction) 
is treated as a quantity of 0{e) leading to extra terms in the nucleon Lagrangian which cancel 
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for a detailed example of such 



the contributions causing this behavior. We refer to Ref. 
a decoupling construction. 

Irnplementing this decoupling of the A resonance near the chiral limit leads to the re- 
sult [3|: 



{9\ 



\Z 



mz 



+ 



167r2F2 



+ 4 <; C^^^(A) + 



155 
972 



91 



17 



36 



9a 



+ 7 



SSE 



In 



X 



277rF2A, 



277r2F2 



l-^ln.^ 



+ 8i§^(25,,-57,°)^n 



2An 



m. 



(9) 



10 
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N 



N S N 



(2) 



(3) 



(4) 




A S iV 




(5) 



(6) 



(7) 
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FIG. 3: Diagrams contributing to the quark-mass dependence oi gA up to 0{e^) in the SSE scheme. 



with 



^,SSE _ ^ ^2 ^ 1 2 ^2 „0 /'„0^3 

^ ~ 167r2F2 81 "^^^ " 2 ~ 9 ~ 



^= ^ + (10) 



In the SSE counterterm combination 

C'^^iX) = B;{Xr^ - 2glBl,{\r^ (11) 

A denotes the scale utilized in dimensional regularization, which can be freely chosen as the 
ChEFT results of course do not depend on the employed regularization scheme or scale. 
In order to make this scale independence explicit, the A dependence of all the couplings 
(controlled by their (3 functions) has of course also to be taken into account (see Table 1 in 
Ref. fl for details). 

In Eqs. OH]) and (ITU]) g\ denotes the axial coupling constant of the nucleon in the chiral 
limit. The parameters and Aq are the pion decay constant (with the physical value 
92.4 MeV) and the real part of the A^A mass splitting, whereas ca and gi denote the leading 
axial A^A and AA couplings, respectively. It is understood that these four parameters are 
also taken in the chiral limit. 

Throughout this work we have converted the quark-mass dependence of gA into a de- 
pendence on the mass of the lowest lying 0~ boson (i.e. the "pion") in the theory via the 
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relation 



2 5om[l + C(m)] 



(12) 



where m = {rriu + m^) /2 denotes the average of the u and d quark mass and the parameter 
Bq = —{0\qq\0)/F^ measures the size of the SU(2) chiral condensate in the chiral hmit. 
The leading term in Eq. (fT2l) corresponds to the GOR- relation 29|. The indicated higher 
order quark-mass terms only start contributing beyond the 0{e^) considered here. We also 
note that Eq. f[T2]) is consistent with the available lattice QCD data for the quark-mass 
dependence of the pion mass (see, e.g., Fig. 1 in Ref. [301]). 
Here we employ Eq. ([9]) in a different representation: 



gA (oo) 



9'a + 



ml , nx9 1 



(9 



\3 

a) 



In 



167r2F2 3247r2F2 



+ 



1087r2F2 



X 



277rAoF2 



+ [25clgiAl - hlc\g\Al - 2Ac\g\ml] 



A| Ao 

— — — arccos 

8l7r2F2Ao 



^ 254^1 (2A2 - m2) 2Ao , c\g\{3ml - 38Al) , 2A 



1627r2F2 ' bAn^F^ + ^(O • (13) 

In Eq. (11 3p we have analytically continued the logarithms of Eq. (jH]) to the region m^^ > Aq, 
because that is the region where our Monte Carlo data lie. Furthermore we have introduced 
the couplings Bg{X) and B2q{X) used in heavy baryon chiral perturbation theory (HBChPT), 
which can be calculated from the SSE couplings via the relations 

1154^1 



B^iXr'^ 



b;{x) 



BUX) 



19447r2F2 
354 



25cl^i ^^2Ao 
6487r2F2 ^ A 
2Ao 



r2 
^A 



In 



(14) 
(15) 



4327r2F2 1447r2F2'" A 
Thus we are able to exploit directly the empirical information on these couplings (see Eq. fl2Up 
below) obtained from a HBChPT calculation of vrA^ — * titiN in Ref. 311 • Another advantage 
of utilizing the couplings Bl{X) and BIq{X) becomes clear when we study the chiral limit 
behavior of Eq. (fT3l) : 



9a{oc) 



9\ 



{9\ 



167r2F2 



ml + A[B;iX)-2g'iB;,iX)] 



mz 



(9 



\3 

a) 



(16) 



In Eq. ( IT6|) it becomes manifest that the couplings defined in Eqs. (fT^ and (fT5|) ensure 
that the 0{e^) SSE result of Eq. (fT3l) displays the same chiral limit behavior as the C(p^) 
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HBChPT result for the quark-mass dependence of qa given in Ref. 32] • In Ref. 3j it was 
made sure that g\ and the leading non-analytic quark-mass dependence ~ In m^r agree in the 
HBChPT and the SSE calculation. Here this mapping of the two field theories is extended 
to the terms ~ m^. A discussion of this point will be given in Ref. 



. B, 3 (i.e. 



However, we want to emphasize again that the 0(e'^) SSE results of Refs 
Eq. (Q) and our Eq. (fT3!l are all equivalent and provide identical results, they only differ in 
the definitions of the employed counterterms. 



IV. CHIRAL EFFECTIVE FIELD THEORY (FINITE VOLUME) 



Applying the methods discussed in Ref. 



14| (for a recent more detailed investigation of 



these finite volume corrections see Ref. 3^) we can now repeat the calculation of qa in 
the SSE scheme for a finite spatial cubic box of length L in order to obtain the volume 
dependence of qa- We define 

gAiL)=gAi(x^)+AgAiL). (17) 

For Qa (oo) we utilize the 0{e^) SSE result of Eq. (fT3l) . From the diagrams of Fig. [3] we 
obtain for AgA (L) to 0{e^) in SSE 

^9a{L) = - — — ^- 





n 
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+ 



i9Af K 
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25c\gi 



cWa 



+ 



8cWa r 



Ko {L\n\f{m^, y)) - ^'"^'^j,"""' K, (L|n|/(m., y)) 



i^o(^N/(m.,y))- ^'^'^j,'"-^^ i^,(L|n|/(m.,y)) 



Ac 



27 -kF^ Ao ^ L\n\m^ 



Kq {L\n\f{m^,y)) 



Ki {L\n\f{m^,y)) 
L\n\f{m„,y) 



(18) 



with 



/(m^, y) = ^ml + + lyA^ . (19) 

The S indicates the omission of the \n\ = contribution in the sum over all vectors n 
with integer components. The contributions of the individual Feynman diagrams shown in 
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Fig. [3] to ^qa are given in Appendix A. The axial coupling of the nucleon in a finite volume 
depends on the very same parameters as the infinite volume result. 

Note that the recent analysis of Ref. [^^j is based on the calculation by Beane and 



Savage [?], which does not employ decoupling constraints in the chiral extrapolation function 
of gA{oo). Presumably this results in a fit function where the coupling g\ does not agree with 
the value of qa at vanishing pion mass. This in turn could lead to different predictions for 
the finite volume corrections, since the analytical formula for AgA{L) should be independent 
of the use of decoupling constraints but is expressed in terms of the same parameters as the 
infinite volume result. In Appendix A we attempt a comparison between our result for 
A(7^(L) of Eq. (ITS!) and the calculation by Beane and Savage 

Equation fll7l) now allows us to access a larger set of lattice data than the formula f[T51) . 
without introducing new parameters. 

V. FIT RESULTS 

As in previous applications of SSE results to the chiral extrapolation of baryon properties 



(see, e.g., Ref. 



28 



35| ) we do not have enough data points at sufficiently small masses 
to fit all parameters. So we must fix some of the parameters at reasonable values. What 
is known from phenomenology? The analysis of (inelastic) vrA^ scattering, in particular the 



p.oce=. ..N M, -ggests that choos.ng the physical pion o,a. a= the =cale A one 

has [3| 

59^(A = mfy^) = (-1.4±1.2)GeV-2 , Bl^{\ = mf^') ^ ^ . (20) 
The coupling B'^q cannot be observed independently of Eg-, as it accompanies a structure in 



the chiral Lagrangian which is proportional to the equation of motion (see e.g. Ref. 36|). 
Hence the separation between the two couplings given by Eq. ( l20i) can be utilized without 
imposing any model assumptions. 

Therefore we set A = 0.14 GeV in the following and take i?2o t)e zero in order to utilize 
this valuable information from scattering experiments. Furthermore, from analyses of tiN 



scattering and vr-photoproduction in the A(1232) resonance region one knows 37| the real 



part of the A^A mass splitting to be 0.271 GeV at the physical point. We note that a recent 



analysis 38|] of the quark-mass dependence of the (real part of the) mass of A(1232) found 



a slightly higher value for the splitting between the mass of the nucleon and its first excited 
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TABLE IV: Fit results. 



fit 


9a 




91 


xVdof 


A 


1.15(12) 


-0.71(18) GeV-2 


2.6(8) 


1.41 


B 


1.26(8) 


-0.62(3) GeV"2 


3.3(7) 


1.54 



state in the chiral limit, Aq ~ 330 MeV, albeit with a large uncertainty. We will study 



the implications of this finding in a future communication [33|]. The imaginary part of the 
complex mass of A (1232) is also known 37|: In an O(e^) analysis within SSE it corresponds 
to a strong decay into ttN intermediate states with the coupling value ca = 1-5 y. In 
the previous analysis of the quark-mass dependence of qa of Ref. a smaller value for 
Ca was utilized. Ultimately such an issue can only be decided by performing simultaneous 
global fits to several observables sensitive to A(1232) contributions within the same ChEFT 
formalism. 

At the physical pion mass we also know the very precise value for the axial coupling 
^exp _ ]^ 2695 of the nucleon from neutron beta decay analyses 371], which in the 0{e^) 



SSE analysis of Ref. together with a set of quenched lattice data led to the estimate 



g'^ ^ 1.2. Little is known about gi. In the SU(6) quark model one would expect gi = 



|l.2 = 2.16. For the pion decay constant ChPT analyses suggest a chiral limit value of 



86.2 MeV 



40|. In the following we fix Aq = 0.271 GeV, ca = 1.5, = 86.2 MeV and 



leave g^, gi, Bg = Bg{\ = 0.14 GeV) as fit parameters. 

To be on the safe side as far as the applicability of ChEFT is concerned we start with a 
fit that takes into account only the six data points at our lowest quark masses (data points 
3, 6, 12, 13, 14 and 17 in Tables HI and HB . Indeed, we cannot expect ChEFT to work at 
pion masses well above 500 — 600 MeV, because the pion-loop integrals start to become more 
and more sensitive to scales beyond A^^, ~ 47?^,^ for such large masses (see the discussion in 



Ref. 



4l| for details). Fortunately, we can include one (/5, k) pair where we have simulations 



for three different volumes. 

Thus we deal with a three-parameter fit of six data points, which we shall call Fit A. 
In the fit function we have, of course, to include the finite volume correction, i.e. we use 
gA^L) = gA{oo) + AgA{L) with gA{oo) given in Eq. (fT3|) and AgA^L) taken from Eq. (ITSl) . 
The resulting values for the fit parameters are given in Table IIVI 
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In Fig. m we plot the data with the finite size correction AgA{L) subtracted together 
with the fit curve. If the fit would describe the volume dependence of the data perfectly, 
the data points from simulations differing only in the volume (simulations 12, 13 and 14) 
would fall on top of each other in the plot. As can be seen in Fig. Hlthis is indeed rather 
well satisfied within the error bars. Furthermore, subtracting /S.gA{L) has moved the "large 
volume" results at our smallest quark masses from all four (3 values even closer together, 
corroborating our previous impression that lattice artefacts are small. 



1.4 
1.2 

< 

1.0 



0.8 

0.0 0.1 0.2 0.3 0.4 0.5 

m,' [GeV'] 

FIG. 4: Fit A of the qa data. The finite size correction has been subtracted from the simulation 
results. Only data points included in the fit are shown. The star represents the physical point. 

Remarkably enough, our fits do not break down when data at somewhat higher pion 
masses are included. Indeed, if we take into account also the data from simulations 9, 10 
and 11, where m^r ~ 0.77 GeV, we obtain the results labeled as Fit B in Table HVl They are 
well compatible with the numbers from Fit A. In Fig. [5] the fit curve is confronted with the 
data. 

Using the fit parameters from Fit B, we plot in Fig. [6] the dependence of qa on the box 
size L for m^r = 0.594 GeV together with our three data points at /? = 5.29, n = 0.1359 and 
for ttItt = 0.769 GeV along with the three data points at 13 = 5.29, k = 0.1355. In addition 
we show the behavior to be expected at m^r = 0.35 GeV and at the physical pion mass. 
These plots show clearly how well Fit B reproduces the volume dependence of our data. 
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FIG. 5: Fit of the data over an extended range of pion masses (Fit B). The finite size correction 
has been subtracted from the data. In contrast to Fig.Hl only the resulting infinite volume numbers 
are shown for the two masses at /3 = 5.29. The three data points around ~ 0.7 GeV^ have not 
been included in the fit. 

In Fig. [7] we display the volume dependence of the data and the fit in yet another way. 
There we plot our results as they were obtained in the respective volumes versus m^. For 
the curves in this plot we take the parameters from Fit B. The finite L values correspond 
to the volumes used in the simulations 12, 13, 14. 

In order to study the influence of the smaller volumes (simulations 9, 10 and 12, 13) 
we have repeated Fit B without these small-volume data, i.e. using only the results from 
simulations 3, 6, 11, 14 and 17. This yields g'^^ = 1.21 ± 0.19, 5^ = (-0.62 ± 0.06)GeV-^ 
and Qi = 2.8 ± 1.7. Of course, the errors have increased, but the fitted parameters are fully 
compatible with the outcome of Fit B. 

In order to estimate the uncertainty caused by the ambiguities inherent in the scale 
setting procedure (see Sect. [TTl) we have not only employed tq = 0.467 fm in the fits but 
also tq = 0.5 fm. Moreover, we have considered the two possibilities of working with ro/a 
evaluated at the respective quark mass and using the chirally extrapolated value of ro/a. 
For we find numbers between 1.12 and 1.21 in Fit A, while they vary between 1.19 and 
1.33 in Fit B. The results for (in GeV"^) lie between -0.60 and -0.83 in Fit A and 
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FIG. 6: Results for (7^ at /? = 5.29, k = 0.1359 (corresponding to a pion mass of about 0.6 GeV, 
upper plot) and at /? = 5.29, k = 0.1355 (corresponding to a pion mass of about 0.77 GeV, lower 
plot) plotted versus L. The full curves have been computed for the corresponding pion masses using 
the parameters from Fit B. The dotted (dashed) curves show the volume dependence expected from 
these parameters for 771,7^ = 0.35 GeV (mjr = 0.14 GeV). 
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FIG. 7: Results for qa in the simulated volumes together with fit curves from Fit B. 



between —0.54 and —0.72 in Fit B. The values for gi range between 2.1 and 3.2 in Fit A 
and between 2.5 and 4.1 in Fit B. We shall take half of these spreads as our estimate of the 
systematic error due to setting the scale. 

Let us now discuss the values that our fits have yielded for the parameters. For the axial 
coupling in the chiral limit g\ we find a result below the value at the physical point. This 
is in agreement with the analyis of Ref. based on the quenched data of QCDSF 



42|. 



The small negative value for the coup 
analyses of inelastic vrA^-scattering 3|, 



ing Bg is also entirely within expectations from the 
3l| (see Eq. [20]) . Finally, we note that the value for 



the axial A A coupling gi in this new analysis is much closer to the SU(6) quark model result 



than the number found in Ref. 



Bl], reassuring us that the contributions of A(1232) to the 



axial coupling of the nucleon are under control. The main cause for this more reasonable 
value of gi is the larger value for the coupling ca utihzed in the present study 43| . 

We note that the curves for gA in a finite volume presented previously in Ref. 10| — which 
were consistently above our infinite-volume curve — can be obtained from the present analyis 
by setting all A(1232) contributions in A(yf^(L) identically to zero. The complete finite 
volume effects found here within 0{e^) SSE, in contrast, are negative at large quark masses, 
consistent with the findings from the simulation. The resulting finite-volume corrections 
therefore arise from a destructive interference between A^vr and Avr loop effects. This point 
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will be discussed in more detail in a future publication 33|]. Such a strong cancellation 



between Goldstone Boson loops around octet and decuplet baryon intermediate states has 



also been observed in Ref. [4^ in an SU(3) HBChPT calculation of the axial current of 
the nucleon in infinite volume, albeit under the additional assumption of degenerate octet- 
decuplet multiplets. 

VI. SUMMARY 

We have evaluated the axial coupling constant of the nucleon in lattice QCD on gauge 
field configurations generated with two degenerate flavors of dynamical non-perturbatively 
improved Wilson fermions. The necessary renormalization of the axial vector current has 
been performed non-perturbatively within the RI'-MOM scheme, except for the mass- 
dependent factor: 6^ has been calculated in tadpole improved boosted perturbation theory. 

For two k) combinations we have performed simulations on three different volumes 
allowing us to study finite size effects. For the chiral extrapolation as well as for the de- 
scription of the volume dependence we have made use of ChEFT. For this purpose we have 
rewritten the expression for the quark-mass dependence of qa derived within the framework 
of the SSE in Ref. in a form which is particularly suitable for our application, and we 
have extended the calculation to cover also the volume dependence. With this formula at 
hand, we could perform a fit to our data obtained on lattices of different spatial extent. 
Note, however, that we had to fix some parameters at phenomenologically reasonable val- 
ues. Then finite volume effects are reproduced surprisingly well down to box lengths of 
about 1 fm. Nevertheless, simulations at smaller quark masses will be necessary to confirm 
our findings^ 

In Ref. 



45| it has been argued that, in the chiral limit, 1/3 of the axial charge of the 
nucleon is to be found at infinite distance from the nucleon, due to a delocalization effect. 
Furthermore, it was suggested that this phenomenon could lead to large finite volume effects. 
While there is agreement on the delocalization phenomenon, it has been disputed that this 
effect will cause calculations of qa in a finite system to have particularly large finite volume 



effects 461]. Notwithstanding, our fits suggest, for example, that at the physical pion mass 
the axial charge is reduced by a factor ^ 2 if the nucleon is confined to a periodic box of 
extent L ^ 1.5 fm. This result is in good agreement with a model calculation [9|, in which 
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the effect is attributed to periodic boundary conditions rather than delocahzation of the 
axial charge. 

Given that the SSE is based upon the long distance pion dynamics around nucleon and 
A(1232) matter states, we conclude that contributions from A(1232) are crucial in under- 
standing the quark-mass (and volume) dependence of qa- This confirms the findings of 
Ref. 3| and the old expectations (at the physical point) based on the observation that the 
integral in the Adler-Weisberger sum rule, which represents g\ — 1, is dominated by the 
A(1232) resonance. The values for the couplings involving A(1232) dynamics employed 
in our chiral extrapolation curve are consistent both with scattering phenomenology and 
with the expectations of the SU(6) quark model. In particular the value for the axial AA 
SSE coupling gi ^ 3.0 in this new analysis is much more consistent with expectations from 
phenomenology than the value obtained in Ref. 

Inclusion of simulation data at several values of the box length L is crucial to increase the 
number of data points at sufficiently small pion masses. The two-dimensional surface (m,r, L) 
allows for a much better determination of the individual values of the effective couplings. 
With the new dynamical simulation data presented here and our 0{e^) SSE analysis we can 
extract the chiral limit value from the extrapolation without including any constraints 
on the physical point. As our final result we quote g^ from Fit B, which takes into account 
pion masses below 800 MeV. We find 

.7° =1.26(8)(7). (21) 

The first error is statistical, while the second error is an estimate of the systematic uncer- 
tainty caused by the ambiguities in setting the scale. Repeating the fit with the value of gA 
at the physical pion mass as fit parameter instead of g\ we obtain 

^7^(mf^^) = 1.31(9)(7). (22) 

Ultimately, one would like to determine all parameters of ChEFT solely from a fit to 
lattice data. Such an enterprise would require a joint fit of results from simulations with 
dynamical quarks for many static nucleon and A(1232) observables, which presently is out 
of reach due to the paucity of lattice data for pion masses below 600 MeV in many of these 
observables. 
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APPENDIX A: AMPLITUDES IN FINITE VOLUME 



The amplitudes of the contributions of the eight Feynman diagrams of Fig. [3] to qa in 
the infinite volume can be found in the appendix of Ref. [2!. Here we give the individual 
contributions of the eight diagrams to Aqa defined in Eq. ( TT71) : 
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Additional details are given in Ref. 43|]. Note that the finite- volume shift in diagram 1 
arises from the Z factor of the nucleon. In Eqs. flAll) - flA6P 5^ denotes the Pauli-Lubanski 
vector of ChEFT, whereas corresponds to an external axial-vector background source 
interacting with the hadronic system; u(ri), n(r2) are the non-relativistic spinors of the 
incoming/outgoing nucleon with 4-momenta r^, respectively. The bilinear combination 
riW^f] with b = 1,2,3 denoting the isovector index of the background source encodes the 
isospin dependence of the current. The sums extend over all vectors n with integer compo- 
nents excluding n = 0. For further details on the notation we refer to Ref. 2|. The couplings 
and parameters occurring in these relations are defined in the main text. 

In TableEwe have atten,pted to .elate ow coupUngs A, c.. ,u A, and K to the o„e= 
used by Beane and Savage in Ref. [7|. However, due to the decoupling constraints of the SSE 
scheme (see the discussion in Sec. IIIip the correspondence between the two sets of couplings 
shown in Table |V] is strictly true only at leading order in the chiral expansion. 

While we are also employing a different representation of the finite-volume shifts in terms 
of the Bessel functions, numerically we can reproduce the results shown in Ref. j3] by utilizing 
the two parameter sets discussed there and shown in Table I VII However, the sign of the 
finite-volume shift observed in the Monte-Carlo simulation of qa reported in this work does 
not agree with either of these two coupling scenarios, which would both lead to an increase 
of qa when the volume is decreased for pion masses above 140 MeV. Finally, we emphasize 
again that our set of couplings of the effective theory in the SSE scheme has been obtained 
for the complete result gA{L) of Eq. ( fTTl) . and not just for the finite- volume shift AgA{L) as 
discussed in Ref. LDl. At O(e^), differences between the SSE scheme used in this work and 

Q 

the approach of Ref. [7|] — aside from numerical differences in the couplings — only manifest 
themselves in the (?yi(oo)-part of gA{L). 
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